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ABSTRACT 

The detection and characterization of the first transiting super-Earth, CoRoT-7b, has required an 
unprecedented effort in terms of telescope time and analysis. Although the star does display a radial 
velocity signal at the period of the planet, this has been difficult to disentangle from the intrinsic 
stellar variability, and pinning down the velocity amplitude has been very challenging. As a result, 
the precise value of the mass of the planet — and even the extent to which it can be considered to 
be confirmed — have been debated in the recent literature, with six mass measurements published so 
far based on the same spectroscopic observations, ranging from about 2 to 8 Earth masses. 

Here we report on an independent validation of the planet discovery, using one of the fundamental 
properties of a transit signal: its achromaticity. We observed four transits of CoRoT-7b at 4.5 /im 
and 8.0 /im with the Infrared Array Camera (IRAC) onboard the Spitzer Space Telescope, in order to 
determine whether the depth of the transit signal in the near-infrared is consistent with that observed 
in the CoRoT bandpass, as expected for a planet. We detected the transit and found an average depth 
of 0.426 ± 0.115 mmag at 4.5 nm, which is in good agreement with the depth of 0.350 ± 0.011 mmag 
(ignoring limb darkening) found by CoRoT. The observations at 8.0 /xm did not yield a significant 
detection. The 4.5 /zm observations place important constraints on the kinds of astrophysical false 
positives that could mimic the signal. Combining this with additional constraints reported earlier, we 
performed an exhaustive exploration of possible blends scenarios for CoRoT-7b using the BLENDER 
technique. We are able to rule out the vast majority of false positives, and the remaining ones are 
found to be much less likely than a true transiting planet. We thus validate CoRoT-7b as a bona- fide 
planet with a very high degree of confidence, independently of any radial-velocity information. Our 
Spitzer observations have additionally allowed us to significantly improve the ephemeris of the planet, 
so that future transits should be recoverable well into the next decade. 

In its warm phase Spitzer is expected to be an essential tool for the validation, along the lines of 
the present analysis, of transiting planet candidates with shallow signals from CoRoT as well as from 
the Kepler Mission, including potentially rocky planets in the habitable zones of their parent stars. 
Subject headings: binaries: eclipsing — planetary systems — stars: individual: CoRoT-7 — stars: 
statistics — techniques: photometric 



1. INTRODUCTION 

Among the known exoplanets, a few special cases 
stand out as the objects that inaugurated the study 
of the physics of Earth- like exoplanets. CoRoT-7b 
(jLeger et al.ll2009t ) is the first super-Earth for which the 
mass and radius have been estimated, and provided the 
first real constraints on models of the formation, struc- 
ture, and evolution of small and potentially rocky exo- 
planets. Kepler-lOb has recently been announced as the 
first rocky planet found by the Kepler Mission (Batalha 
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et al. 2010), with a mass of Mp = 4:.6~^\''^ and a radius 

of Rp = 1.4:16~^q'^IqRq. We now also have examples of 
interesting planets that are intermediate in mass and ra- 
dius between the Earth and Neptune, such as GJ 1214 b, 
with Mp = 6.55 ± 0.98 Mm (Charbonneau et al. 2009) 
and Rp = 2.64 ± 0.13 i?e (|Berta et al.ll2011[ ). 

The first detection of a super-Eart h was made p ossible 
by the successful CoRoT satellite (jBaglin et all [20021. 
The discovery was made in the course of observations 
in the first long run of this mission in the direction of 
the anti-center of the Galaxy (LRaOl), which took place 
from 2007 October to 2008 March. A smaU transit-like 
signal was identified with a depth of 0.35 mmag and du- 
ration of 1.3 hr, recurring with a period of 0.8535 days 
and being consistent with a super-Earth size planet or- 
biting a bright (V = 11.7, K = 9.8) G9 dwarf star. The 
discovery triggered a series of follow-up observations to 
clarify the origin of the shallow transit signal. The de- 
tection an d the follow- up campaign have been fully de- 
scribed bv lLeger et all |2009). Bruntt et al. (2010) sub- 
sequently performed a detailed spectroscopic analysis of 
the CoRoT-7 star and determined an improved stellar 
radius of i?* = 0.82 ± 0.04 Rq. This resulted in a revised 



2 



p lanet radius o f i^p = 1.58 ± 0.10 i?®. 

iQueloz et al.l (|2009D reported on an extensive obser- 
vational campaign carried out with the HARPS instru- 
ment on the ESO 3.6 m telescope at La Silla, with the 
goal of detecting the Doppler signal of this small ob- 
ject and measuring its mass. CoRoT-7 is an active star, 
however, with starspot-induced photometric variability 
at the '^2% percent level modulated by the stellar ro- 
tation period of about 23 days. The radial velocity of 
the star is dominated by an irregular signal with an am- 
plitude several times larger than the sought-after sig- 
nature of the transiting planet. Under these circum- 
stances, measuring the mass accurately is a challenging 
task that depends strongly on how the activity-induced 
variability is handled, on assumptions about possible ad- 
ditional non-transiting planets that may be contribut- 
ing (up to two have been cons idered), and on the atti- 
tude toward systematic errors. IQueloz et ahl |2009 ) pro- 
duced the first mass estimate of Mp = 4.8 ± 0.8 M®. 
Subsequent authors have reported different values us- 
ing the same HARPS observations or subsets thereof. 
These estimates are not always consistent with each 
other within their formal errors, and v ary consider- 
ably in significance le vel: 6.9 ib 1.4 M m (jHatzes et al.l 
\20m. 2.3 ± 1.5 Ma. (IFont et al.l [20lTh. 8 ± 1.2 Mm 
(iFerraz-Mel lo et al.l W ill). 5.7 ± 2.5 Ma, _ (jBoisse et alj 



I2011D . and 7.4±1.2M^ (jHatzes et al.ll201lD . The impact 
of the discrepancies is not insignificant, as these estimates 
lead to rather different mean densities and therefore dif- 
ferent internal structures for the planet as inferred from 
current theory. While some of the more recent deter- 
minations appear to favor a higher mass for the planet, 
at the 2(7 level the estimates range all the way from 
to 10 Mffi , and for some of the lower values (iPont et al.l 
120111 : iBoisse ct al. 2011) the statistical significance of the 
Doppler detection of CoRoT-7b is considerably less com- 
pelling. 

Similar difficulties are expected to occur for other can- 
didates with shallow transits, in which "confirmation" 
in the usual sense of the word by the detection of a 
radial- velocity signature that is at the limit of detectabil- 
ity with current instrumentation will be very challenging. 
With the recent announcement of a large number of shal- 
low transit candidat es discovered by the Kepler satellite 
(jBorucki et al.|[20Tl|) , obtaining assurance that these sig- 
nals correspond to bona-fide super-Earth size planets, as 
opposed to a false positive, is among the most urgent 
tasks that lie ahead, and is of primary importance for 
the statistical interpretation of the results. 

Motivated by the lingering problems with the CoRoT- 
7 b mass determination described above, and the impli- 
cations for the robustness of the detection, we use this 
case here to illustrate the application of a powerful tech- 
nique to help "validate'O shallow transit candidates in- 
dependently of any radial-velocity information. It makes 
use of near-infrared observations with the Spitzer Space 
Telescope, and is based on the premise that true tran- 
sits are achromatic signals, to first order, so that the 

In the context of this paper "confirmation" as used above 
refers to the unambiguous detection of the gravitational influence 
of the planet on its host star (e.g., the Doppler signal) to establish 
the planetary nature of the candidate; when this is not possible, 
we speak of "validation" , which involves an estimate of the false 
alarm probability. 



transit depth as observed in the near-infrared should 
be the same as in the CoRoT passband (ignoring the 
effects of limb darkening). If the candidate is the re- 
sult of a blend, however, the depth can be significantly 
different. As an example, if the observed transit were 
due to an object eclipsing a background star of 0.5 Mq 
and sun-like age and metallicity, then its transit in the 
Spitzer 4.5 ^m bandpass would be 3.2 times deeper than 
in the CoRoT bandpass. With its infrared passband, 
Spitzer affords the maximal wavelength separation from 
the CoRoT photometry (passband around 650 nm), and 
places very strong constraints on possible false positive 
scenarios, as we describe below. 

Even with these constraints, and others available from 
follow-up observations carried out and reported by the 
CoRoT team, it is not possible to rule out all possi- 
ble ble nd configura t ions for CoRoT-7b, as recognized 
also bv ILeger et al.l ( |2009D . Thus, the main goal of this 
work is to more exhaustively explore the wide variety of 
false positive scenarios that can mimic the light curve, 
to obtain a realistic estimate of the blend frequency that 
may be expected. We aim to provide an independent 
assessment of the confidence level that the signal is of 
planetary nature. Rather than focusing solely on the 
likelihood of a blend (frequentist approach) , as in previ- 
ous studies, we adopt a Bayesian approach in which we 
compare the blend frequency with a prior for the like- 
lihood of a planet (odds ratio). To evaluate the blend 
frequency we mak e use of the BLENDER technique in- 
troduced by Torres et al.l (12004 [2011). wi th further de- 
velopments as described bv iFressIn etall (l20Til) . This 
methodology has been applied successfully to validate 
a number of shallow transit signals from the Kepler 
Mission includin g Kepler-9 d, Kepler- 10 c, Kepler- llg, 
and KeDler-19 b (ITorres et al. 20 11: Frcssin et al. 201 1|; 
' Lissauer et al.ll2011l : iBallard eraI.1f2Clllh . 

We begin by describing our Spitzer observations (Sec- 
tion [5]), and then briefly summarizing the data used here 
along with other follow-up observations relevant to this 
investigation (Section[3]). This is followed by the BLENDER 
analysis that examines the vast space of parameters for 
false positives by synthesizing realistic blend light curves 
and comparing them with the CoRoT photometry (Sec- 
tion H]). We next estimate the expected frequency of 
blends and compare it with the expected frequency of 
planets (Section [5]). As shown below, this analysis is 
able to validate CoRoT-7b as a planet without relying 
on any radial- velocity information. 

2. Spitzer OBSERVATIONS OF CoRoT-7B 

The Infrared Array Camera (IRAC; 'Fazio et al."l200^ 
of the Spitzer Space Telescope (Werner et al. 200'^~ob- 
tains simultaneous images in four bandpasses. A 5.'2 x 
5'.2 field of view (FOV) is imaged in one pair of band- 
passes (3.6 and 5.8 ^m), and a nearly adjacent FOV 
imaged in the second pair (4.5 and 8.0 /xm). The two 
blue channels employ InSb detectors, whereas the red 
channels use Si:As IBC detectors. All four arrays are 
256 X 256 pixels. While the present Warm Spitzer mis- 
sion is restricted to the two shorter wavelengths, the 
data discussed in this section were obtained just prior 
to the spacecraft entering that phase in May of 2009. 
We elected to monitor CoRoT-7 in only one channel pair 
(4.5 and 8.0 /itm), as even if the stellar flux in these band- 
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passes is slightly smaller, we have obtain ed more precise 
observations in the past using this pair (|Knutson et alj 
l2008t iFressin et al.ll200"9l) . It is also farther in wavelength 
from the CoRoT bandpass, and the expected depth dif- 
ference if the signal comes from a blend scales with the 
wavelength difference, as we describe later. 

We used IRAC to observe the primary eclipse of 
CoRoT-7b on UT 2009 April 22, 23, 24 and 25, obtain- 
ing data at 4.5 and 8.0 /xm. We were able to observe it 
in full array mode in both channels for a total duration 
of 17.0 hours, including four transits lasting 75 minutes 
each. We observed the target in the IRAC stellar mode, 
in which the camera gathers simultaneously two 10.4 s 
integrations at 4.5 and 8.0 fj,m. Therefore, we gathered 
respectively 1075, 1212, 1212, and 1212 images both at 
4.5 and 8.0 /im during the four transits we observed. Our 
goal was to detect the transit signal at 4.5 fim, but as the 
8.0 fim observations are simultaneous and automatic, we 
present them hereafter for completeness. We describe 
below our observations in two sections, as the InSb de- 
tectors used for IRAC channels at 4.5 fim require a dif- 
ferent treatment than the Si:As detectors of the IRAC 
8.0 fim channel. 

We were mindful that the signal we were looking for 
was at the limit of what was possible to obtain with a 
few transits from Spitzer, as it has mainly been used to 
look at brighter stars and deeper eclipses. Our expecta- 
tion in terms of statistical significance for the detection 
of a single transit was 2.1cr at 4.5 /im (4.2(t when com- 
bining the four transits), scaled on the detection level we 
achieved in our previous studies of TrES-1, TrES-3, and 
TrES-4 (iCharbonneau eTal] 120081: iKnutson et ahl 120081: 
IFressin et al.ll2009D . 

2.1. J^.d fim observations 

The contribution of the background to the total flux 
from CoRoT-7 is low in the 4.5 fim IRAC bandpass, con- 
tributing only 0.35% to the total flux in an aperture with 
a 5-pixel diameter centered on the position of the star. 
We obtained the lowest RMS time series using an aper- 
ture with a radius of 5.0 pixels. We allowed the position 
of our aperture to shift with the position of the star in 
each image. We estimated the background in each image 
from an annulus with an inner radius of 12 pixels and an 
outer radius of 20 pixels centered on the position of the 
star. 

We determined the position of the star in each im- 
age by fitting a two-dimensional Gaussian to the position 
of the star. Agol et al. (2010) have recently completed 
a comparative analysis of different methods to estimate 
the stellar centroid in Spitzer images. Their best results 
are obtained by fitting a two-dimensional Gaussian to 
the star PSF, compared to a fiux-weighted centroiding, 
and parabolic fitting. This Gaussian algorithm uses the 
7x7 pixel sub-array from the image centered on the 
brightest pixel. It fits a two-dimensional Gaussian to 
this array, allowing its center, amplitude, and width to 
vary. It then uses a non-linear Levenbcrg-Marquardt al- 
gorithm to optimize these parameters (Mar kwardt 2009,). 
We compared this technique with the position-weighted 
sum of the fiux in a 4-pixel radius disk centered on the ap- 
proximate position of the star, which we previously used 
on similar magnitude targets ("Charbonn eau et al.l 120081 : 
lKnutson"et al...2008; .Frcssin et al..,2009, ). The Gaussian 



fit proved to be better regarding two criteria. First, 
the scatter between the position estimates was 1.3 times 
smaller. Second, the scatter of the differential position 
between our target star and a nearby reference star (2 
magnitudes fainter in channel 2 and at 28") was smaller 
and did not show any correlation with the intra-pixel 
position. 

The dominant instrumental systematic effect in the 
first two IRAC b andpasses is due to a well-known intra- 
pixel sensitivity fReach et al.' '2005"; Charbonncau et "ah 
2005^^ 2008; Moralcs-Caldcron ct al. 2006; Knutson et al. 
20081 ). Fluxes at these two wavelengths show a strong 
correlation with the intra-pixel position of the star on the 
detector, at a level comparable to the expected depth of 
the eclipse. We used the following parameters to fit the 
observed flux as a linear function of the subpixel position: 

/ = /o * [ci + C2{x - xo) + C3{y - yo)] , (1) 

where fo is the original flux from the star, / is the mea- 
sured flux, X and y denote the location of the Gaussian-fit 
centroid of the star on the array, xq and yo sue the coor- 
dinates of the center of the pixel containing the peak of 
the star's point spread function, and C1-C3 are free pa- 
rameters in the fit. We excluded the in-transit measure- 
ments, based on the known ephemeris, in order to avoid 
suppressing the transit depth. For the x and y positions 
above we calculated the centroid of the target in each im- 
age and found that the pointing jitter was around 0.12 
pixels (Of 14) over the course of a visit in both x and y. 
The pointing drift of the telescope appears to occur on a 
longer timescale than the exposures. Therefore, instead 
of using the actual position estimate at individual expo- 
sure times, we smoothed the x and y curves as a function 
of time and used the smoothed position instead. 

In contrast t o previous observations of HD 189733 
and H P 209458 ()Knutson et al. 2008: Charbonncau et all 
l2008f l. we found that adding quadratic or higher-order 
terms to this equation, or even cross-terms, did not im- 
prove the fit significantly, likely due to the lower signal- 
to-noise ratio of the present observations. 

After correcting for the intra-pixel sensitivity, a de- 
creasing trend was still visible that is likely to be an in- 
strumental effect related to the detector or telescope, and 
has also been seen in observations of TrE S-3 and TrES-4 
(|Knutson et al.ll2008t fFressin et aLl l2009^. two stars with 
similar brightness. We corrected for this effect by fit- 
ting the data in both channels with a linear function 
of time. This term was fitted simultaneously with the 
correction for the intra-pixel sensitivity, so that we can 
accurately characterize the additional uncertainty in the 
depth and timing of the eclipse introduced by these cor- 
rections. That is to say, we solved for four parameters 
including a constant term, a linear function of x posi- 
tion, a linear function of y position, and a linear func- 
tion of time. We also trimmed the first 30 minutes of 
data that show a larger scatter, as we have done in pre- 
vious cases. The fit was performed with a Markov Chain 
Monte Carlo (MCMC) method (Ford 2005; Winn et al.l 
I2007D with 5 x 10^ steps, where we set the uncertainty of 
the individual measurements equal to the standard de- 
viation of the out-of-transit data after correction for the 
various detector effects. 

Prior to the fit we carried out an initial trimming 
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Fig. 1. — Transit light curves of CoRoT-7b in channel 2 (4.5 /xm) 
observed on UT 2009 April 22, 23, 24 and 25, with best-fit curves 
representing instrumental effects overplotted in red. Data have 
been binned in 7.3 minute intervals, then offset by a constant for 
the purposes of this plot. 



within our aperture, discarding outliers farther than 3.5cr 
from the local median flux (defined as the median of a 
15-minute window centered on the data point). We also 
removed measurements for which the identified position 
of the photocenter x or y deviated by more than 3.5ct 
from the same 15-minute median position. This global 
trimming excludes 6.6% of the data points in the four vis- 
its in the 4.5 /im bandpass. We excluded outliers greater 
than 3.5(7 during each step of the chain, as determined 
using the residuals from the model light curve, from our 
evaluation of the function. We rescaled the value of 
the function to account for the fact that we are vary- 
ing the number of pixels included in the fit. 

After producing the chain, we searched for the point 
in the chain where the value first falls below the me- 
dian of all the values (i.e., where the code had first 
found the best- fit solution), and discarded all the steps 
up to that point. Figured] shows the four individual light 
curves and the respective fits for the instrumental effects 
that we removed from these curves before looking for the 
transit signal itself. 

Next we carried out a second Markov chain fit of the 
transit signal itself on the trimmed data. We initially al- 
lowed the individual depths and times of the four transits 
to vary independently, along with the normalized semi- 
major axis a/Ri,, and the inclination angle of the orbit. 
We calculated o u r tran sit curves using the formulation by 
iMandel fc Agoll ()2002[ ). AUhough smaU at 4.5 fim, limb 
darkening was taken into account i n our modeli ng ac- 
cording to the four-parameter law jjy I Claret' ('2000|), with 
coefficients taken from the work of i Sing (2010) and stellar 
parameters Tcff = 5250 ± 6 K, logq = 4 47 ± .05, and 
[Fe/H] = -HO. 12 ± 0.06 from lBruntteTall (poToh . While 
the individual depths and transit time parameters were 
consistent with transits occurring at the expected period 
from the CoRoT light curve, the signal-to- noise ratios 
were too poor to provide a meaningful constraint on the 
geometric parameters. We therefore chose to restrict the 



TABLE 1 

Best-fit transit parameters 



Parameter 


CoRoT^' 


Spitzer 4.5 ^tm 


P (day) 


0.853585 ± 0.000024 


0.853590 ± 0.000006 




4.27 ±0.20 




Inclination (deg) 


80.1 ± 0.3 




Depth (mmag) 


0.350 ±0.011 


0.426 ±0'. 115 



" ILeger et all l[2009h . 



space of free parameters for the Markov chain analysis 
to the transit depth, a/i?^,, inclination angle, and period, 
an d we held the trans it epoch fixed at the value reported 
bv lLeger et all (|2009D . 

Non-random sources of noise in transit and eclipse pho- 
tometry — such as instrumental systematics and stellar 
variability — could dominate the error budget in the de- 
rived system parameters. This is true for ground-based 
data, and also turns out to be true for space-based data. 
The higher stability of space measurements is offset by 
the fact that smaller effects are being measured, and cor- 
respondingly smaller level s of random error are achieved 
by collecting more signal. iPont et all (|2006l) showed how 
neglecting this type of noise could lead to an underes- 
timate of the actual uncertainties by a large factor. In 
the present study, a realistic error estimate is important 
and we therefore attempted to assess the possible im- 
pact of non-random noise. The presence of correlated 
noise is obvious in our raw photometric sequence. We 
corrected for it to first order, but obviously the correc- 
tion cannot be perfect. We used the single-parameter 
description of the correlated component of the noise pro- 
posed by [Po nt et al. (2006), with the further simplifica- 
tion of I Winn et al.l (|2009 l adapted to regularly-sampled 
data, to estimate the impact of the residual systemat- 
ics after decorrelation. We repeated the MCMC anal- 
ysis by using modified uncertainties a^^^ = cr^ + ncr^, 
where cr^ and are the 'white' and 'red' components 
of the noise (i.e., random and correlated, respectively), 
and n is the number of data points during the eclipse. 
We estimated a^- from the dispersion of the flux after 
decorrelation obtained with slightly different, reasonable 
decorrelation procedures. We used ar = 0.000108 mag as 
the dispersion between three different decorrelation tech- 
niques (i.e., fitting the intra-pixel answer and the transit 
simultaneously, adding quadratic terms to the intra-pixel 
answer, and the decorrelation technique previously de- 
scribed). 

Table [T] collects our results, and the Spitzer time series 
is shown in Figure [2] along with the fitted model. The pa- 
rameters we derive for the planet are in good agreement 
with those based on the light curve as observed by the 
CoRoT satellite. In particular, the near-infrared and op- 
tical depths are consistent within the errors. Also shown 
in Figure [2] is a model based on the optical lightcurve 
parameters, for comparison with the Spit zer curve. Al - 
though we do not reach the precision Le ger et al.l (poM ) 
obtained for the geometric parameters based on their de- 
tection of 153 individual transits, we were able to improve 
the precision in the period determination of the planet 
significantly, as our Spitzer observations were gathered 
some 640 planetary orbits after the original CoRoT data. 
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Fig. 2.— Transit light curve of CoRoT-7b in the IRAC 4.5 /xm 
bandpass after removal of instrumental effects, folded on the ex- 
pected period. Data have b een n ormalized to remove detector ef- 
fects (see discussion in Sect. \2A\i . and binned in 9.2 minute inter- 
vals. The black curve is the best transit model fit to the Spitzer 
light curve, imposing the CoRoT ephemeris. The transit depth 
is 0.426 it 0.115 mmag. The green curve shows the model ex- 
pected for the supe r-Earth tr a nsit s cenario, adopting the CoRoT 
parameters from iLeger et al.l II2009I ), which lead to a depth of 
0.350 it 0.011 mmag. The good agreement in the depth indicates 
the transit is achromatic, as expected for a planet. 

To put our results on a more quantitative basis, we 
computed the reduced values corresponding to the 
case where no transit at all is shown in the data, and 
the case of a transit with a depth corresponding to a 
planet such as CoRoT-7b. We find that the super-Earth 
scenario gives a reduced (1.0039) very close to the best 
fit to the data (1.0033), whereas the no-transit scenario 
provides a poor fit (1.0276). 

Since the transit depth we measure in the near in- 
frared is about the same as in the optical, this argues 
against blends composed of stars of much later spectral 
type. Importantly, we can use the error bar on the mea- 
sured depth to rule out blends involving stars of a dif- 
ferent temperature that would have produced a signifi- 
cantly different depth in infrared. Indeed, the depth of 
a blend in the Spitzer bandpass compared to the optical 
depends solely on the infrared contribution of the blend- 
ing star to the total flux relative to its contribution in 
the CoRoT bandpass. For example, if the relative con- 
tribution in the infrared of a background eclipsing binary 
is three times larger in infrared than in the visible, then 
the observed transit in the infrared would be three times 
deeper. The relative infrared-to-visible contribution of 
a background star is also independent of its distance to 
the target, to first order (i.e., ignoring the effects of inter- 
stellar dust). As a consequence, for a given target star, 
the observed depth increase or decrease of a blend when 
observed in the Spitzer bandpass is only a function of 
the spectral type of the contaminating star. Using an 
isochrone of solar age and metallicity representative of 
the background stars in the Galactic plane, it is straight- 
forward to compute the change in the relative infrared-to- 
visible flux contribution of such a star as a function of its 
mass, and hence the depth increase or decrease in the 4.5 
and 8 /im bandpasses that a blend would display. With 
the t arget propert i es (ma ss M — 0.91 Mq, age « 2 Gyr) 
from iBruntt et all (|2010t ). we find for CoRoT-7 that all 
blended stars below 0.69 Mq would show a transit deeper 
than the 3-a upper limit observed with Spitzer at 4.5 /im 
(0.77 mmag). We make use of this important constraint 
later in Sect. [5] to eliminate a large fraction of potential 



false positives for CoRoT-7b. 

2.2. 8.0 fim observations 

For the observations at 8.0 /xm we used the "preflash" 
technique (Knutson et al. 2009), in which we pointed the 
telescope for 30 minutes towards a bright H II region be- 
fore observing CoRoT-7. This was completed in order to 
reduce the amplitude of the detector "ramp" at 8.0 /im, 
effectively pre-loading the pixels on which the target star 
would be pointed. 

Previous secondary eclipse studies (e.g., Knutson et al. 
2008) have shown that PSF-fitting can provide a better 
signal-to-noise ratio at longer wavelengths. At longer 
wavelengths the flux from the star is smaller and the zo- 
diacal background is larger; we find that the background 
contributes 20% of the total flux in a 3-pixel aperture at 
8.0 ^m. 

At 8.0 /im, we found that the relative scatter in the 
time series after model fitting from the PSF fits was 15% 
higher than in the time series from aperture photometry 
with a 3.0 pixel radius. As a result of this increased scat- 
ter, which is likely produced by discrepancies between the 
model PSF and the observed PSF, we concluded that 
aperture photometry is also preferable. We compared 
the time series using apertures ranging from 3.0 to 4.5 
pixels and found consistent results in all cases, but with a 
scatter that increases with the radius of the photometric 
aperture. 

Previous observations (e.g.. iKnutson et aI]|2009D have 
shown that the ramp is well described as following an 
asymptotic shape, with a steeper rise in the first 30 min- 
utes of observations. We corrected for this effect by fit- 
ting our time series with the following function; 

/ = /o*[ci+C2ln(dt)] , (2) 

where fo is the original flux from the star, / is the mea- 
sured flux, and dt is the elapsed time in days since the 
start of the observations. 

As the expected statistical significance of the signal at 
8.0 /im was very small (i.e., 1.1 cr), we chose to fix the 
transit time and period of the transit to values found in 
our 4.5 /im study and verified that the residual signal in 
the four visits at 8.0 /xm was compatible with the tran- 
siting planet scenario. We carried out a Markov Chain 
Monte Carlo fit to the data as described in Sect. 12.11 si- 
multaneously fitting Eq. [5] and a transit model, with the 
depth as the only free parameter in the latter. No signif- 
icant correlations were found between the variables. As 
a further check we repeated these fits adding a quadratic 
term of the form In(dt) in Eq. [21 and found that the value 
of for our best-fit solution was similar to the previous 
value, so we chose not to include this additional term. 

Figure [3] shows the detrended (from the ramp), 
trimmed, binned, and folded channel 4 light curve. We 
estimated the best-fit transit depth using a Markov 
Chain Monte Carlo fit in the same way as we did for 
the 4.5 /im data, and found an eclipse depth of 0.11 ± 
0.30 mmag. The noise level is too large to say anything 
about a dip corresponding to the planet eclipse, but it 
disallows a blend scenario involving a significantly red- 
der star. Specifically, proceeding in the same way as 
described at the end of Sect. I2.1[ a star of 0.5 Mq would 
be excluded at the 3a level. 
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Fig. 3.— Transit light curve of CoRoT-7b in the IRAC 8.0 Atm 
bandpass after removal of instrumental effects, phased with the ex- 
pected period and with the best-fit transit curve overplotted (see 
text). Data have been normali zed t o remove detector effects (see 
discussion in Sect. [2T1 and Sect. [2T2t . and binned in 9.2 minute in- 
tervals. The black curve represents the best transit model fit adopt- 
ing the CoRoT ephemeris, which gives a statistically insignificant 
depth of 0.11 it 0.30 mmag. The green curve shows the fit expected 
for the sup e r-Ear th transit scenario, using the parameters from 
ILeger et al.1 1 |2009|'| . 

3. CoRoT PHOTOMETRY AND FOLLOW-UP 
OBSERVATIONS 

The photometric observations used in the next sec- 
tion to investigate blends are essentially the same as 
used by ILeger et al.l (|2009[ ). with a somewhat different 
de-trending of the data. We started with the N2 (sci- 
ence grade) light curve, discarding the initial few days, 
which had a time sampling of 512s, and worked only 
with the r emainder, where the time-sampling was 32s 
(see ILeger et al. 2009L for a detailed description of the 
N2 light curve). We identified and clipped outliers using 
a running 5-point median filter, and then modeled the 
intrin sic stellar variabihty u sing an iterative non-linear 
filter (|Aigrain fc IrwinI 120041 ) . This consists in applying 
a 5-point boxcar filter followed by a median filter with a 
width of 1 day, 3cr clipping the residuals, and iterating 
until no more points are clipped. The resulting slowly- 
varying component was then subtracted from the original 
light curve, to give the time-series used in this work. We 
made use of the white-light data only. Of the extensive 
follow-up observ ations carrie d out by the CoRoT team 
as described by ILeger et "aLl (|2009[) , the most relevant 
for our study are the high-resolution imaging observa- 
tions obtained with the NACO instrument on the VLT. 
These data exclude any blends capable of mimicking the 
signal outside of 0'.'4 from the target down to 6.5 mag 
fainter than the target. Additionally, the spectroscopic 
observations reported by these authors allowed them to 
rule out companions that would be bright enough to be 
visible in the spectrum (as a second set of lines). In par- 
ticular, their near-infrared (i^-band) spectroscopy with 
CRIRES on the VLT rules out most companions brighter 
than about 7% of the flux of the target in that band 
(Am « 2.9). For the analysis below we adopt a more 
conservative limit of Am = 1.0. 

4. BLENDER ANALYSIS 

Blend scenarios satisfying the CoRoT observations 
were explored with the BLEND ER technique (|Torres et al.l 
l200ll20TllFressin et al.ll20lll) . by examining the quality 
of the fit to the CoRoT-7 photometry of a large array of 
synthetic model light curves. These synthetic light curves 



result from the combined light of three objects compos- 
ing the blend: the main target, and an unresolved eclips- 
ing pair along the line of sight with the same period as 
detected in CoRoT-7, whose eclipses are attenuated by 
the light from the target. The eclipsing system may be 
physically associated (hierarchical triple), or may be in 
the background or foreground, and the eclipsing pair may 
consist of two stars, or a star transited by a larger planet. 
BLENDER uses the detailed shape of the transit light curve 
to weed out scenarios that lead to the wrong shape for 
a transit. The properties of the three obj ects were taken 
from model isochrones bv lGirardi et al.l ([20001 in order 
to synthesize light curves, and these artificial light curves 
were compared with the observations in a sense (i.e., 
by computing the sum of the squared residuals normal- 
ized by the photometric uncertainties). For practical rea- 
sons blend light curves were calculated here in the Kepler 
passband, which is sufficiently similar to the passband of 
the CoRoT satellite for our purposes, in both the cen- 
tral wavelength and width[3 The properties of the main 
star (essentially its intrinsic brightness, which affects the 
dilution of the eclipses) were constrained by the spec- 
troscopic observations of CoRoT-7 (Tcff = 5250 ± 60 K 
log q = 4.47 ± 0.05 [Fe/H l = -t0.12 ± 0.06: IOueloz et all 
l200a IBruntt et al1l2010l ). and held fixed. Those of the 
other two objects (referred to here as the 'secondary' for 
the eclipsed star and 'tertiary' for the eclipsing object) 
were varied over a wide range and parametrized by mass. 
For hierarchical triples the stellar components were con- 
strained to lie on the same isochrone, while for back- 
ground/foreground blends the properties of the intruding 
star were taken from a representative solar-metallicity, 
3 Gyr isochrone. 

The distance between the binary (dEe) and the 
target (dcoRoT-y) was expressed for convenience in 
terms of the distance modulus difference {AS — 
5 log [dEB/rfcoRoT-r]), which is independent of interstel- 
lar extinction, and varied between —5 and +10. The im- 
pact parameter was allowed to vary between and unity, 
and the mass of the secondary stars ranged from 0.1 Mq 
to 1.4 Mq. For stellar tertiaries we explored masses from 
0.1 Mq up to the mass of the secondary; for tertiaries 
that are planets (assumed to be dark) we allowed their 
sizes to be up to 2.0i?jup. We restricted our blend sim- 
ulations to circular orbits for all eclipsing pairs, as no 
stellar or planetary systems with periods as short as that 
of CoRoT-7b are known to have eccentric orbits, nor are 
expec ted to from theoretical arguments (see, e.g.. lMazehl 
l2008f ). Synthetic light curves were generated with a de- 
tailed eclipsing binary code, including proximity effects 
(tidal and rotational distortions), limb-darkening, grav- 
ity brightening, and contamination from third light. Dif- 
ferential extinction was included for chance alignments. 

Blends providing poor fits compared to the CoRoT-7 
photometry in a sense were considered to be ruled 
out. This enables us to place constraints on the kinds 
of objects composing the eclipsing pair that yield viable 
blends (i.e., acceptable fits), including their size or mass, 
as well as other properties of the blend such as the over- 
all brightness and color. For further details and applica- 
tions of BLENDER to other transitin g planet can d idates , 
we refer the reader to the work of iTorres et al.l (|2011[) . 

(see |http: //keplergo ■ arc ■ nasa.gov/CallbratlonResponse ■ shtml[ l 



7 



iLissauer et"all (|2011[ ). and lFressin et alj ()201lD . 



4.1. Background or foreground stars transited by a 

planet 

We consider first the case of blends involving a back;- 
ground or foreground star falling within the CoRoT aper- 
ture, and transited by a larger planet. Because of the 
very short 0.8535-day orbital period of CoRoT-7 evolved 
stars (giants or subgiants) are ruled out, so we focus in 
the following on main sequence stars. We find that there 
is a very large range of spectral types (masses) permitted 
for the secondary star, shown in top panel of Figure 01 
as well as a wide range of relative distances between the 
eclipsing pair and the target. This is indicated by the 
darker areas in the figure, delimited by the white con- 
tour representing light curve fits that differ from the best 
transiting planet model by a difference corresponding 
to a 3-tT confidence level. However, other constraints 
available for CoRoT-7 strongly restrict the number of 
these false positives. This is illustrated in the bottom 
panel of the figure, in which these additional constraints 
are superposed on the same blend landscape displayed 
in the top panel. In particular, by comparing the pre- 
dicted r — Kg color of each blend against the measured 
color of the star from Exodat {r - Ks ^ 1.723 ± 0.025; 
iLeger et al.l [20091 ) . we find that a large fraction of the 
secondaries with significantly different masses than the 
primary are ruled out because the blends would be too 
red or too blue compared to the known color index of 
CoRoT-7 (by more than 3a). These excluded regions 
are indicated by the blue hatched areas in the bottom 
panel of Figured Other blends are excluded because the 
secondary star would be very bright (within one magni- 
tude of the target), and would have been noticed spec- 
troscopically, if unresolved in the high-resolution imag- 
ing described earlier. The section of parameter space 
excluded by this brightness criterion is indicated by the 
green hatched area, limited from above by the green di- 
agonal line corresponding to a magnitude difference of 
1.0 mag. Note that, as mentioned in Sect. [3l this is a 
very conservative brightness limit as a large fraction of 
stars within 2 or even 3 magnitudes of CoRoT-7 would 
most likely have been detected spectroscopically. As indi- 
cated before, our Spitzer observations exclude secondary 
stars less massive than approximately 0.69 A/q, which 
are common constituents of background blends. This 
additional constraint is shown by the shaded gray area 
leftward of 0.69 M0. There is some overlap between the 
gray region and the blue and green hatched areas, indi- 
cating that those blends are excluded by more than one 
observational constraint. 

An example of a blend that provides an acceptable 
fit to the CoRoT photometry is shown in Figure [7] (top 
panel). This scenario (location marked with a cross in 
Figure mi is not excluded either by its color or its bright- 
ness, but is clearly ruled out by our Spitzer observations 
because it produces a much deeper transit in the near- 
infrared (see bottom panel of Figure [7]), which would 
have been easily detected. To summarize, the combi- 
nation of the Spitzer, color, and brightness constraints 
removes many but not all blends involving a background 
or foreground star transited by a larger planet. Those 
that remain reside in the area of Figure |4] (bottom) la- 
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Fig. 4. — Top: Map of the surface (goodness of fit) for 
blends involving background or foreground stars transited by a 
larger planet. The vertical axis represents the distance between 
the background pair of objects and the primary star, cast for con- 
venience in terms of the difference in the distance modulus AS (note 
that A<5 is not equivalent to the magnitude difference between the 
background system and the main star because of the effects of dif- 
ferential extinction, which are included in our simulations). Only 
blends inside the solid white contour match the CoRoT light curve 
within acceptable limits {3cr, where (j is the significance level of the 
difference; see lFressin et al"]|2011l) . Lighter-colored areas (red, 
orange, yellow) mark regions of parameter space giving increasingly 
worse fits to the data (4(7, 5(T, etc.), and correspond to blends we 
consider to be ruled out. Bottom: Same diagram as above with 
the addition of observational constraints from follow-up measure- 
ments. The constraint from our Spitzer observations is represented 
by the shaded gray area to the left of 0.69 A/q; all false positives 
with secondary masses smaller than 0.69 Mq can be rejected, as 
they lead to transit depths at 4.5 /im that are inconsistent with our 
measurements. Blends in the hatched green area are also ruled out 
because they are bright enough to be detected spectroscopically 
(Am < 1.0 mag, represented by the solid green line). The hatched 
blue regions correspond to blends that can be excluded as well be- 
cause of their r — Kg colors, which are either too red (left) or too 
blue (right) compared to the measured value for CoRoT-7 hy more 
than 3cr (0.075 mag). The combination of all of these constraints 
leaves only a reduced area of parameter space (labeled 'Allowed 
Region') where blend models give tolerable fits to the CoRoT light 
curve, and are not ruled out by any of the follow-up observations. 
These blends are all brighter than Am = 5.5 mag (dashed green 
line). The white cross to the left of the Allowed Region marks the 
location of a representative blend ruled out by Spitzer, for which 
the predicted light curves are shown in Fig. [7] (see text). 



beled "Allowed Region" . 

4.2. Background eclipsing binaries 

For the case of a background eclipsing binary com- 
posed of two stars, interestingly we find that no combi- 
nation of relative distance and stellar properties for the 
eclipsing pair gives an acceptable fit to the CoRoT light 
curve. The reason for this is that all such blend con- 
figurations that can potentially reproduce the detailed 
shape of the transit also lead to out-of-eclipse brightness 
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Fig. 5. — Similar to Fig. |4] (and with the same color scheme) 
for blends involving background or foreground stars transited by 
another star. Using BLENDER we find that no such background 
eclipsing binary is able to reproduce the observed photometry at 
a level better than 8.3 a (white contour) in comparison to the 
of the best planet fit. Therefore, this excludes all background or 
foreground eclipsing binaries consisting of two stars as potential 
false positives. 

changes (elhpsoidal variations) with an amphtude that 
is not seen in the data, and that are a consequence of 
the very short orbital period. This imphes that back- 
ground blends of this kind can be confidently ruled out, 
an imp ortant conclusion t hat does not follow from the 
original TLeger et al.l ()2009[ ) analysis. Figure [S] shows that 
the blends that yield the best fits to the photometry are 
excluded at the 8.3 cr level or higher. This result is sig- 
nificant, as it substantially reduces the overall likelihood 
of blends for CoRoT-7b. 

We note that white dwarfs are also excluded as poten- 
tial tertiaries. Although the range of their radii over- 
laps with those of the planets considered earlier in this 
section, their considerably larger mass would once again 
induce significant ellipsoidal variation in the light curve, 
which is ruled out by the observations. As was the case in 
the previous section, giant stars do not constitute viable 
secondaries because of the very short 0.8535-day period 
of the signal. 

4.3. Hierarchical triples 

Finally, for eclipsing binaries (consisting of two stars) 
that are physically associated with the target in a hier- 
archical triple configuration, we find that the blend light 
curves invariably have the wrong shape to mimic a true 
transiting planet signal, for any combination of stellar 
parameters for the secondary and tertiary. Either the 
depth, duration, or the steepness of the ingress/egress 
phases of the transits provide a very poor match to the 
CoRoT photometry, resulting in differences compared 
to the best transit model corresponding to several hun- 
dred a. These scenarios are therefore all excluded. On 
the other hand, if we allow the tertiaries to be planets, 
the blend fits are somewhat better over a wide range of 
masses for the secondaries when transited by a planet of 
the appropriate size, but not quite at the 3cr level or bet- 
ter. The best for ^ blend of this type corresponds to 
about a 3.4cr departure from a planet model. Although 
our formal 3-cr limit is reasonable and consistent with 
common practice, it is still somewhat arbitrary and one 
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Fig. 6. — Similar to Fig. [4] for the case of hierarchical triple 
systems in which the secondary is transited by a planet. After 
taking into account the Spitzer constraint (gray shaded region on 
the left), as well as the constraints on the r — Ks color and those on 
the brightness of the secondaries from spectroscopy (blue and green 
hatched regions, respectively), we find that all triple configurations 
are excluded as false positives. 



may argue that fits that are only marginally worse might 
still be tolerable. Even accepting this possibility. Fig- 
ure [5] shows that all of these blends are excluded by a 
combination of constraints from Spitzer, color index, and 
brightness, even those that are lOtr or more away from 
the quality of a planet model. 

4.4. Summary of blends 

From the simulations described above, the only vi- 
able blend scenarios for CoRoT-7b are those involving a 
larger planet transiting a background or foreground star. 
BLENDER restricts those blends to the area labeled "Al- 
lowed Region" in Figure 01 These configurations involve 
stars between 1.0 and 5.5 mag fainter than the target 
within the CoRoT aperture. 

5. VALIDATING CoRoT-7B 

The a priori frequency of stars in the background or 
foreground of the target that are orbited by a transit- 
ing planet and are capable of mimicking the photomet- 
ric signal may be estimated from the density of stars in 
the vicinity of CoRoT-7, and the frequency of transiting 
planets with the appropriate characteristics. The rele- 
vant area around the target is that in which stars of each 
brightness would go undetected in the high-resolution 
imaging reported by the CoRoT team. To obtain the 
number density (stars per square degree) we in ake use 
of the Galactic structure models of iRobin et al.l ([2003), 
and we perform this calculation in half-magnitude bins, 
as illustrated in Table [5] For each bin we further re- 
strict the star counts using the constraints on the mass 
of the secondaries supplied by BLENDER (see Figure S]). 
These mass ranges are listed in column 3, and the result- 
ing densities appear in column 4(3 Bins with no entries 
correspond to brightness ranges excluded by BLENDER. 

We note that the precipitous drop in the numbers listed in 
column 4 below a magnitude of 16.7 is not due to a real decrease 
in stellar density (which in fact rises for fainter stars), but to the 
fact that the mass range allowed for blends in this magnitude bin 
is significantly reduced. 
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For the maximum angular separation (pmax, column 5) 
at which stars of each brightness would escape detection 
we have ado pted t he value 0'.'4, based on the report by 
iLeger et alJ (|2009l ) that companions outside this range 
would have been seen in their VLT/NACO image down 
to a magnitude difference of 6.5 mag compared to the tar- 
get. We note that this Pmax value is a conservative limit, 
as stars at closer separations would have been detected 
if they had smaller magnitude differences, which would 
reduce the blend frequency. The result for the number 
of background stars in each magnitude bin is given in 
column 6, in units of 10~^. 

To estimate the frequency of transiting planets that 
might be expected to orbit these stars ( and lead to a 
false positive) we rely on the results from iBorucki et all 
([2011), who reported a total of 1235 planet candidates 
among the 156,453 Kepler targets observed during the 
first four months of the Mission. These signals have not 
yet been confirmed to be caused by planets, and there- 
fore remain candidates until they can be thoroughly fol- 
lowed up. However, the rate of false positives in this 
sample is expected to b e quite small (10% or less; see 
IMorton fc" Johnson I 2011D . so our results will not be sig- 
nificantly affected by the assumption that all of the can- 
di dates are plane t s. We further assume that the census 
of IBorucki et al.l ()201lD is largely complete. After ac- 
counting for the additional BLENDER constraint on the 
range of planet sizes for blends of this kind (tcrtiaries of 
0.24-1.42 i?jup), we find that the transiting planet fre- 
quency is /planet = 571/156,453 — 0.0036. Multiplying 
this frequency by the star counts in column 6 of Table [H 
we arrive at a total blend frequency (BF) listed at the 
bottom of column 7, BF = 4.2 x 10"'^. 

This figure represents the a priori likelihood of a false 
positive, and we note that it is approximately 3 orders of 
magnitude smaller than indicated by the calculations of 
ILeger et"al] (|2009D . However, we do not consider this to 
represent the "false alarm rate", as the expected likeli- 
hood of a transiting planet of the characteristics implied 
by the transit signal is also very small. We adopt here 
a Bayesian approach analogous to that employed to vali- 
date previous Kepler candidates, in which our confidence 
in the planetary nature of the signal will depend on how 
the blend likelihood compares to the a priori likelihood of 
a true transiting planet (PF), addressed below. Thus, we 
seek to estimate the odds ratio PF/BF. This is a signif- 
icant conce ptual difference com pared to the frequentist 
approach bv lLeger et al.l ()2009() . who considered only the 
likelihood of a blend (BF). 

Implicit in the blend frequency calculation above 
(BF = 4.2 X 10~^) is the 3a criterion on the quality of 
the light curve fit relative to a transit model fit that we 
used as a condition for a blend scenario to be acceptable 
(see Sect.|4]). For a fair comparison, we use a similar 3cr 
criterion to establish the a priori transiting planet fre- 
quency (PF), for the numerator of our odds ratio. We 
estimate PF by cou nting the Kepler candidates in the 
IBorucki et al.l ()201lf ) sample that have radii within 3cr 
of the value determined from t he best fit to the Co RoT- 
7 data (Rp = 1.58 ± 0.10 i?®; iBruntt et alllMol) . We 
find 231 candidates within this range, giving a planet 
frequency PF = 231/156,453 = 0.0015. 

Thus, the likelihood of a planet is more than 3,500 
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Fig. 7. — Example of a blend scenario (solid curves) that re- 
produces the visible light curve from CoRoT (black crosses), but 
that is ruled out solely by our Spitzer observations (and not by 
any of the other follow-up observations). It involves an MOV star 
of 0.53 -Mq with V = 17.0 in the background of the target star, 
eclipsed by a Jovian planet with a radius 1.4 times that of Jupiter. 
Top: Folded light curve, full phase. Bottom: Enlargement around 
the transit of CoR,oT-7 b. The dashed curve shows that the transit 
depth predicted for this same blend scenario in the near-infrared 
4.5 fim bandpass is much greater than in t he o ptical, which is in- 
consistent with our findings (see also Sect. [2?T] and Fig.[2)l. 

times greater than that of a false positive (PF/BF — 
0.0015/4.2 X 10""^), which we consider sufficient to in- 
dependently validate CoRoT-7b as a true planet with 
a high degree of confidence. We note that our blend 
frequency calculation assumes the 1,235 candidates cat- 
aloged bv IBorucki et al.l (|2011l ) are all true planets. If we 
were to assume conservatively that as many as 50% are 
false positives in the radius range specific to the tran- 
siting planet case (an unlikely propositi on that is also 
inconsistent with other evidence; sec How ard et alll2010l: 
IBorucki et al.ll20lll) . the planet likelihood would still be 
nearly 1,800 times greater than the likelihood of a blend, 
implying a false alarm rate sufficiently small to validate 
the signal. 

Additionally, in estimating the a priori true planet fre- 
quency (PF) and the frequency of larger planets involved 
in blends, we have not placed any restriction on the peri- 
ods of either kind of planet. One may expect, for exam- 
ple, that limiting the periods of larger planets involved in 
blends to be similar to that of CoRoT-7b, which is quite 
short, might reduce the blend frequency quite consider- 
ably. On the other hand, a similar period limit imposed 
on the true planet frequency would decrease PF as well, 
though these effects may not be equal. Thus, it is pos- 
sible that the odds ratio would be altered. As it turns 
out, our calculations without any period constraints are 
the most conservative. For example, limiti ng the periods 
to be shorter than 5 days, the census of IBorucki et al.l 
()2011|) indicates that larger planets are comparatively 
less common than the smaller ones, resulting in an im- 
provement in the PF/BF odds ratio of a factor of 2.2. 
Restricting the periods to be less than 3 days, the odds 
ratio improves by a factor of 2.8. 

6. DISCUSSION 

Our observation of four transits of CoRoT-7b with the 
Spitzer Space Telescope at a wavelength of 4.5 /itm has 
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resulted in the detection of the transit with a depth of 
0.426±0.115nimag, which is con sistent with t he Co RoT- 
7 b planet scenario described in lLeger et all ()2009D . Al- 
though the signal-to-noise ratio of this infrared detec- 
tion is relatively low, it has been obtained with a now 
standard treatment of the Spitzer data with no a pri- 
ori knowledge of the transit parameters (aside from the 
timing windows in order to plan the four Spitzer visits). 

Warm Spitzer is currently the only facility available 
that has the capability of detecting such shallow tran- 
sits at wavelengths that are sufficiently separated from 
the CoRoT (or Kepler) passbands to be helpful. In this 
case the observations were successful, and the transit at 
4.5 fim is shown to have virtually the same depth as in 
the optical. This places a strong constraint on the color 
of potential blends, which are restricted to have secon- 
daries of similar spectral type as the primary star. In 
the case of CoRoT-7b this allows us to rule out most 
cool stars (which constitute the majority of background 
stars) as potential contaminants. 

The detailed analysis of the CoRoT-7 photometry with 
BLENDER combined with constraints from other obser- 
vations eliminates the vast majority of possible blend 
scenarios. This includes all background eclipsing bina- 
ries composed of two stars, most of the scenarios in- 
volving chance alignments with a star transited by a 
larger planet, and all hierarchical triple configurations. 
The remaining scenarios are much less likely (by a fac- 
tor of 3,500) than a true transiting planet, thereby val- 
idating the planetary nature of the signal with very 
high confidence. We point out that this conclusion has 
been reached with very conservative assumptions regard- 
ing some of the observational constraints. For exam- 
ple, we have ignored th e fact that t he hi gh-resolution 
VLT/NACO imaging bv iLeger et"all (120091 permits the 
detection of background stars closer than our adopted 
limit of Pmax = 0'.'4 from the target if they are rela- 
tively bright. Additional imaging by those authors with 
FASTCAM on the 2.5 m NOT telescope provides even 
tighter constraints that we have not used, reaching down 
to sensitivities of 4 magnitudes fainter than the target 
at O'.'IS. We have also not considered the full poten- 
tial of spectroscopy to rule out closer companions; we 
have assumed only that stars within 1.0 magnitude of 
the target would have been identified, whereas in reality 
the sensitivity of those observations (not only the near- 
infrared CRIRES spectrum, but also those from HARPS) 
is probably much greater (Am « 2-3 mag). Furthermore, 
we have not made use of other information reported by 
ILeger et"aLl (|2009[ ). such as the constraints from the red, 
green, and blue passbands into which the light from the 
CoRoT instrument can be split, which can further limit 
the pool of potential false positives. Incorporating all of 
these constraints can only reduce the blend frequency, 
resulting in an even greater confidence level for the vali- 
dation of CoRoT-7b as the first super-Earth. 

Our Spitzer observations have also significantly im- 
proved our knowledge of the ephemeris of the first tran- 
siting super-Earth. Transits of CoRoT-7b have not been 
observed with any other photometric facility since the 
end of the LRaOl field observations by the CoRoT satel- 
lite. Given the formal uncertainty in the orbital period 



reported by ILeger et"al1 (|2009D , the accumulated error 
in the predicted times of transit up to the date of our 
Spitzer observations is approximately 20 min, or about a 
quarter of the transit duration. Here we have improved 
the period determination by about a factor of four, to 
P = 0.853590 ± 0.000006 days, ensuring the transits wiU 
be recoverable for at least another decade (provided there 
are no physical mechanisms operating to change the pe- 
riod). 

The Kepler satellite has recently released a large 
number of very pr omising shallow transit candidates 
(jBorucki et al.ll201lD . From the high quality of this pho- 
tometry the detailed transit shapes are often very well 
defined, increasing the power of tools such as BLENDER 
that make use of that information to rule out blend 
scenarios that result in poor fits to the data. The 
exquisite relative astrometric precision delivered by the 
Kepler instrument provides important additional con- 
straints from an analysis of the motion of the photo- 
centers of the images in and out of transit (see, e.g., 
iBatalha et al.l |2010). which can typically exclude 90% 
or more of unresolved stars falling within the photo- 
metric aperture that can potentially contaminate the 
flux measurements. Spectroscopy and high-resolution 
imaging contribute further valuable constraints. These 
tools have already been used to validate small signals 
for which radial-velocity confirmation is currently out 
of reach, including possibly rocky planets. Example s 



include the super-Earth Kepler-9 d iTorres et al. '2011J^ 
and the Neptune-s ize planets Kepler-l Oc ( Fressin et al 
l20ll. Kepler-11 g (iLissauer et aLlbofTl ). and Kepler-19b 
( Ballard et ani20m 

Spitzer observations provide a very effective way of rul- 
ing out a large fraction of the blend scenarios involving 
a background star with a different effective temperature 
(or color) than the target. These chance alignments are 
typically the most serious concern regarding the nature 
of the transit signals. Among the most interesting candi- 
dates expected from the Kepler Mission are rocky planets 
in the habitable zone of their parent stars. The typically 
longer periods of these objects mean that methods of 
confirmation relying on the dynamical influence of the 
planet and/or the quality of the phase- folded photome- 
try (obtained by summing individual transit events) will 
be more problematic, as their efficiency scales down with 
orbital period. We anticipate that the Spitzer telescope 
in its warm phase will prove critical for validating such 
objects, as its efficiency does not depend on period, but 
relies instead on a different intrinsic property of the tran- 
sits, which is their achromaticity. 
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TABLE 2 

Blend frequency estimate for CoRoT-7 b. 









Blends Involving Planetary Tertiaries 


Magnitude Range 


Am 


Stellar 


Stellar Density 


pmax 


Stars 


Transiting Planets 


(mag) 


(mag) 


Mass Range 


(per sq. deg) 


( ) 


(x 10~^) 


0.24-1.42 _Rjup, /planet = 0.36% 






(Mq) 








(xlO-^) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


11.7-12.2 


0.5 












12.2-12.7 


1.0 












12.7-13.2 


1.5 


0.84-1.36 


y 1 


U.4U 


o.Oo 


U.Ulo 


13.2-13.7 


2.0 


0.82-1.33 


137 


0.40 


5.31 


0.019 


13.7-14.2 


2.5 


0.77-1.30 


207 


0.40 


8.03 


0.029 


14.2-14.7 


3.0 


0.69-1.26 


263 


0.40 


10.2 


0.037 


14.7-15.2 


3.5 


0.69-1.22 


383 


0.40 


14.9 


0.054 


15.2-15.7 


4.0 


0.69-1.17 


515 


0.40 


20.0 


0.072 


15.7-16.2 


4.5 


0.69-1.08 


569 


0.40 


22.0 


0.081 


16.2-16.7 


5.0 


0.69-0.99 


637 


0.40 


24.7 


0.090 


16.7-17.2 


5.5 


0.69-0.80 


186 


0.40 


7.21 


0.026 


17.2-17.7 


6.0 












17.7-18.2 


6.5 












18.2-18.7 


7.0 












18.7-19.2 


7.5 












19.2-19.7 


8.0 












19.7-20.2 


8.5 












Totals 






2988 




115.9 


0.42 



Total frequency (BF) = 4.2 X lO"'^ 



Note. — Magnitude bins with no entries correspond to briglitness ranges in wtiicti BLENDER excludes all blends. 
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